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Abstract. Data from the IceCube detector in its 
22-string configuration (IC22) were used to directly 
measure the atmospheric energy spectrum near the 
horizon. After passage through more than 10 km 
of ice, muon bundles from air showers are reduced 
to single muons, whose energy can be estimated 
from the total number of photons registered in the 
detector. The energy distribution obtained in this way 
is sensitive to the cosmic ray composition around 
the knee and is complementary to measurements by 
air shower arrays. The method described extends 
the physics potential of neutrino telescopes and can 
easily be applied in similar detectors. Presented is 
the result from the analysis of one month of IC22 
data. The entire event sample will be unblinded once 
systematic detector effects are fully understood. 

Keywords: atmospheric muons, CR composition, 
neutrino detector 
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Fig. 1: Muon surface energy spectrum measurements 
compared to theoretical models Q. 



I. Introduction 

While the primary goal of IceCube is the detection 
of astrophysical neutrinos, it also provides unique op- 
portunities for cosmic-ray physics (T). One of the most 
important is the direct measurement of the atmospheric 
muon energy spectrum. 

As shown in figure Q] the energy spectrum of muons 
produced in cosmic -ray induced air showers has so far 
been measured only up to an energy of about 70 TeV 
12). The best agreement with theoretical models was 
found by the LVD detector, with the highest data point 
located at = 40 TeV All these measurements 
have been performed using underground detectors. Their 
sensitivity was limited by the relatively small effective 
volume compared to neutrino telescopes. 

With a planned instrumented volume of one cubic 
kilometer, IceCube will be able to register a substantial 
amount of events even at very high energies, where the 
flux becomes very low. The limitation in measuring the 
muon spectrum is given by its high granularity, and 
consequent inability to resolve individual muons. Most 
air showers containing high energy muons will consist 
of bundles with hundreds or even thousands of tracks. 
Since the energy loss per unit length can be described by 
the equation dE jdx = a + bE, low-energy muons will 
contribute disproportionately to the total calorimetric de- 
tector response, which depends strongly on the energy of 
the primary, disfavoring the measurement of individual 
muon energies. 



This problem can be resolved by taking advantage of 
the fact that low energy muons are attenuated by energy 
losses during passage through the ice. In this analysis, 
the emphasis was therefore set on horizontal events, 
where only the most energetic muons are still able to 
penetrate the surrounding material. The primary cosmic 
ray interaction in this region takes place at a higher 
altitude, and therefore in thinner air. The reinteraction 
probability for light mesons (pions and kaons) is smaller 
and the flux of muons originating in their decays is 
maximized. 

The main possibilities for physics investigations using 
the muon energy spectrum are: 

• Forward production of light mesons at high ener- 
gies. While muon neutrinos at TeV energies mostly 
come from the process K — > + X, for kine- 
matical reasons muons originate predominantly in 
pion decays it — > + p, |]4). An estimate of 
the pion production cross section from accelerator 
experiments gives an uncertainty of 

8(a Nn ) ~ 15% + 12.2% ■ log w (E n /500 GeV) 

at x iab > 0.1 above 500 GeV Q. This value 
should also apply in good approximation to the 
conventional (non-prompt) muon flux. 

• Prompt flux from charm meson decay in air show- 
ers [6|. Because of their short decay length, the 
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reinteraction probability for heavy quark hadrons 
is negligible. The resulting muon energy spectrum 
follows the primary energy spectrum with a power 
law index of 7 as —2.7 and is almost constant over 
all zenith angles. Since the non-prompt muon flux 
from lighter mesons is higher near the horizon, this 
means that the relative contribution from charm is 
lowest, and very challenging to detect. 
• Variations of the muon energy spectrum due to 
changes of the CR composition around the knee. 
Since the ratio of median parent cosmic ray and 
muon energy is < 10 at energies J7) above 1 TeV, 
a steepening of the energy-per-nucleon spectrum of 
cosmic rays at a few PeV will have a measurable ef- 
fect on the atmospheric muon spectrum at energies 
of hundreds of TeV. Comparison of the measured 
muon spectrum to various phenomenological com- 
position models was the main focus of this analysis. 

An additional benefit in the case of neutrino detectors 
is that a direct measurement of the muon flux will 
have important implications for neutrino analyses. By 
reducing the systematic uncertainties on atmospheric 
lepton production beyond 100 TeV, the detection po- 
tential for diffuse astrophysical fluxes will be enhanced. 
Also, atmospheric muons serve as a "test beam" that 
allows calibration of the detector response to high- 
energy tracks. 

II. Cosmic Ray Composition Models 

Starting from the hypothesis that most cosmic rays 
originate from Fermi acceleration in supernova shock 
fronts within our galaxy, the change in the energy 
spectrum can be explained by leaking of high energy 
particles. Since the gyromagnetic radius 
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depends on the charge Z of the particle, for a given 
energy nuclei of heavier elements are less likely to 
escape the galactic magnetic field than lighter ones. 

The general expression for the flux of primary nuclei 
of charge Z and energy Eg is 
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where the transition energy E trans corresponds to 
E p ■ Z, E p ■ A or simply E p for rigidity-dependent, 
mass-dependent and constant composition models. The 
parameter e c determines the smoothness of the transition, 
and A7 the change in the power law index. 

Three alternative composition models have been pro- 
posed, which all can be fit reasonabkly well to the total 
cosmic ray flux in the region of the knee [8|. These are: 

• Rigidity-Dependent A7: This is the default com- 
position used in the IceCube downgoing muon 
simulation. It is also the one favored by current 
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Fig. 2: Atmospheric muon energy spectrum at surface 
level averaged over the whole sky as simulated with 
CORSIKA/SIBYLL. 



models of cosmic ray production and propagation 
in the galactic amgnetic field. 

« Mass-Dependent A7: An alternative model that 
also leads to a composition change around the knee. 
The change in the power law index does not depend 
on the charge, but on the mass of the nucleus. The 
best fit proposed in the original paper leads to a 
smaller value for the transition energy and a steeper 
spectrum after the cutoff. 

• Constant Composition: Here, the composition of 
the primary cosmic ray flux does not change. The 
knee is explained by a common steepening in the 
energy spectrum for all primaries occurring at the 
same energy. 

The best measurement of the composition so far was 
done by KASCADE j£|. Its result was consistent with 
a steepening of the spectrum of light elements, but 
depended strongly on the hadronic interaction model 
used to simulate the air showers (SIBYLL or QGSJET). 

The influence of the three composition models on the 
muon energy spectrum is shown in figure [2] While the 
spectrum for the constant composition model gradually 
changes from E~ 3J to E~ 4 , the other two show a 
marked steepening corresponding to the cutoff in the 
energy per primary nucleon. By accurately measuring 
the muon energy spectrum, it is therefore possible to 
significantly constrain the range of allowed cosmic ray 
composition models in the knee region. 

III. Analysis 

The data set used in this analysis is based on the 
IceCube online muon filter, designed to contain all 
track-like events originating from the region below 70°. 
It covers the period from June 2006 to March 2007 
with an integrated livetime of 275.6 days, during which 
IceCube was taking data with 22 strings (IC22). A 
number of quality cuts were applied in order to eliminate 
background from misreconstructed tracks and to reduce 
the median error in the zenith angle measurement to 
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Fig. 3: Relation between energy proxy e sur t and true 
surface energy of most energetic muon in shower. Here 
and in figure |4]the rigidity-dependent composition model 
was used. 



Fig. 4: Simulated muon multiplicity for atmospheric 
showers at closest approach to the center of the Inlce 
detector for different values of e. 



w 0.7°. The final sample corresponded to an event rate 
of 0.146 Hz. 
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TABLE I: Threshold energy for muons passing through 
ice. The energy values correspond to an attenuation of 
99.9%. 



To measure the single muon energy spectrum, it is 
necessary to reduce the background of high-multiplicity 
bundles, whose total energy depends primarily on the 
primary cosmic ray [10|. Since there is no possibility 
to accurately estimate the multiplicity of a downgoing 
muon bundle, the only way to obtain single muons is by 
selecting a region close to the horizon to which muons 
of lower energies cannot penetrate. 

The minimum energy required for muons passing 
through a distance d of ice can be approximated by the 
equation 

E cut (d) = (e bd - l)a/b 

where a = O.lGSGeVm" 1 and b = 0.192 • lO^m.- 1 
ifTTI . The resulting threshold energies corresponding to 
vertical tracks and for tracks at the top and bottom of 
the detector for angles near the horizon are shown in 
Table U 

Two factors determine the upper energy bound of 
this analysis. One is the contribution from atmospheric 
neutrinos, which will eventually dominate the event 
sample at large depths. The other, and more important, is 
the finite zenith angle resolution. Using simulated data, it 
was determined that it effectively limits the measurement 
of the slant depth to a values below 15 km. 




Fig. 5: True surface energy of most energetic muon in 
shower using rigidity-depending composition model for 
different values of e sur f, with fits to Gaussian function. 
All individual distributions are normalized to unity. 



Using the slant depth alone, the range of this analysis 
is therefore insufficient to probe the region beyond 100 
TeV. However, the reach can be extended by incorpo- 
rating information about the energy of the muon as it 
passes through the instrumented volume. 

For muon tracks in the detector, the energy resolution 
approaches Alog 10 (-E) « 0.3 above 10 TeV lfl2l . This 
information can be combined with the slant depth to 
obtain a better estimate for the muon energy at the 
surface. 

A natural way to do this is by defining a surface 
energy proxy e that behaves as 

exp(e sur/ ) oc logn 7 • d s i 

an I 

where n 1 represents the total number of photons 
measured by the detector. Figure [3] shows the re- 
sulting parameter, which has been linerly rescaled in 
such a way that its value corresponds to the mean 
\og{E^ sur f /GeV) for any given bin, provided that 
the muon energy spectrum is reasonably close to the 
standard E~~ 3 7 . 

An important criterion for the applicability of the 
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Fig. 6: Data from one month of IC22 at final cut level 
compared to simulated event rates and fit of empirical 
function exp(a + be + ce 2 ) to constant composition 
distribution. 
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Fig. 7: Ratio of experimental to simulated e sur t distri- 
butions for one month of IC22 data. Uncertainties are 
statistical only and exclude systematic detector effects. 



energy proxy parameter is that over the entire range of 
measurement the muon multiplicity remains low, and the 
influence of high-multiplicity bundles small. Figure [4] 
confirms that this is indeed the case. It should be noted 
here that the most energetic muon typically accounts 
for the dominant contribution to the total energy in the 
detector, such that other tracks in the bundle can be 
neglected. 

The spread in muon surface energies for a given value 
of e sur f is shown in Figure [5] Around the peak the 
distributions can be approximated by a Gaussian whose 
width lies in the range of Alog 10 E ks 0.3 — 0.4. 

IV. Result 

Figure [6] shows simulated event rates in dependence 
of e sur f compared to data at final cut level. Almost over 
the entire range all three models can be approximated 
by the same empirical fit function. Only in the highest 
bin can a distinction be made. 

The experimental data agrees remarkably well with 
the simulation, as can be seen more clearly in Figure 
|7] Despite the steeply falling distribution, the ratio of 



data to simulation remains very close to one over almost 
the entire range. For e > 5, corresponding to E^ > 
100 TeV, the measurent is based on only three data 
events. 

Using the entire year of IC22 data, the predicted event 
yield for 5.1 < e sur / < 5.2 based on the constant 
composition model corresponds to about 10 events. It 
is therefore unlikely, even neglecting systematic uncer- 
tainties, that any of the three models under considera- 
tion could definitely be excluded yet. This situation is 
expected to change as soon as 40-string data can be 
included in the analysis. 

V. Conclusion 

This result demonstrates the potential for an accurate 
measurement of the muon energy spectrum with large 
neutrino detectors. So far only one month of data has 
been considered in the analysis, corresponding to about 
10% of the entire event sample. Nevertheless, the mea- 
surement already covers an energy range almost a factor 
of three above that of the previous upper limit, with very 
good agreement between data and simulation. 

While it will be difficult to make a definitive statement 
about the cosmic ray composition around the knee based 
on IC22 data, it will be possible to confirm the validity 
of cosmic ray air shower models up to previously 
inaccessible energy ranges. 

At the time of writing, the instrumented volume 
of the detector has increased by a almost factor of 
three. Further enlargements are scheduled for the next 
few years. Future measurements of the muon energy 
spectrum will benefit from a larger effective area, and a 
substantial improvement in the angular resolution related 
to the longer lever arm for horizontal muon tracks within 
the detector. 

Once residual systematic detector uncertainties are 
resolved, a comprehensive analysis that accounts for 
both the energy spectrum of individual muons and the 
total shower energy in the detector will be feasible. The 
potential for such a combined measurement is unique to 
large volume detectors. 
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